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Classical simulations of noisy stabilizer circuits are often used to estimate the threshold of a quantum error- 
correcting code. Physical noise sources are efficiently approximated by random insertions of Pauli operators. 
For a single qubit, more accurate approximations that still allow for efficient simulation can be obtained by 
including Clifford operators and Pauli operators conditional on measurement in the noise model. We examine 
the feasibility of employing these expanded error approximations to obtain better threshold estimates. We 
calculate the level-1 pseudothreshold for the Steane [[7,1,3]] code for amplitude damping and dephasing along 
a non-Clifford axis. The expanded channels estimate the actual channel action more accurately than the Pauli 
channels before error correction. However, after error correction, the Pauli twirling approximation yields very 
accurate estimates of the performance of quantum error-correcting protocols in the presence of the actual noise 
channel. 


I. INTRODUCTION 

The threshold theorem of quantum error correction 
promises the accurate implementation of arbitrary size quan¬ 
tum algorithms if the underlying physical errors are below cer¬ 
tain values. The error thresholds depend strongly on the spe¬ 
cific quantum error correcting code, how errors are detected 
and fixed HJ-E1. and what errors are assumed gHZl- Most 
codes have been designed to fix random Pauli errors and error 
correction procedures can be simulated efficiently using the 
stabilizer formalism ED. A broader class of errors including 
Clifford operations Ql and Pauli measurements mo can also 
be included in this formalism. For a single qubit, this extended 
error set has been shown to yield improved approximations of 
realistic error models including amplitude damping ED- 

Here we examine whether these improved approximations 
also lead to more accurate threshold estimates. Specifi¬ 
cally, we calculate the level-1 pseudothreshold for the Steane 
[[7,1,3]] code 02 for two non-stabilizer errors, amplitude 
damping and a depolarization channel along a magic-state 
axis, and compare the exact solution to approximations based 
on Pauli errors or Clifford and Pauli measurement errors. The 
Steane code has been well studied theoretically 1T3UT9I and a 
logically encoded state has been recently demonstrated exper¬ 
imentally li20l . The code is small enough to allow for exact 
simulation similar to recent work on distance-3 surface codes, 
which compared a realistic error model corresponding to Ti 
(amplitude damping) and T 2 (dephasing) processes and an ap¬ 
proximate Pauli error model based on twirling (6J. 

In addition to the pseudothreshold, we are interested in two 
other qualities of the approximation, the accuracy and the hon¬ 
esty. The accuracy is a measure of how close is the state gen¬ 
erated by the approximate evolution to the state generated by 
the exact evolution. We describe an approximation as honest if 
the final state after the approximate evolution is further from 
the initial state than the final state after the exact evolution. 
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In other words, an approximation is honest if it upper-bounds 
the error of the exact evolution. As pointed out by Puzzuoli 
et al. the composition of honest approximations is not nec¬ 
essarily honest ETI . This implies that an approximation that 
is honest at the 1-qubit physical level might lead to a dishon¬ 
est representation of the overall error produced on the system. 
If our goal is to employ our approximate channels to infer 
the performance of error-correcting strategies under realistic 
non-stabilizer noise, then we need to be cautious and be sure 
that they compose in an honest fashion. We provide numeri¬ 
cal evidence that, in the context of an error-correcting circuit, 
an honest approximation at the physical level remains honest 
at the logical level. Furthermore, we show that, for the error 
models studied, physically dishonest approximations based on 
the Pauli channel might lead to approximations at the logical 
level that are both approximately honest and very accurate, in 
agreement with similar results obtained by Geller and Zhou 
ll22l . This suggests that it might not be necessary for the ap¬ 
proximations to be honest at the physical level. 

The paper is organized as follows. In Section [II] we de¬ 
scribe the target realistic error channels and our method for 
generating approximate channels ifTTl . In Section [III] we re¬ 
view the important concepts of honesty and accuracy of an 
approximate channel. In Section|iy]we explain our procedure 
for calculating the pseudothreshold. In Section[V] we present 
our results before concluding in Section [VT] 


II. ERROR CHANNELS 

We review all the error channels introduced in ED- We 
start with the stabilizer expansions to the Pauli channel (PC) 
that are used as models to approximate realistic non-stabilizer 
error channels. Next we discuss two important cases of these 
error channels that lie outside the stabilizer formalism. Finally 
we review two different constraints under which the approx¬ 
imations are performed. All the error channels introduced in 
this section will be expressed in the operator-sum representa¬ 
tion. 
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Throughout the paper we use X, Y, and Z to represent the 
Pauli matrices with associated eigenvectors {|+), |—)}, {| + 
i), | — *)}, and {|0), 11)} respectively. 


A. Efficiently siinulable processes 


One of the main ideas introduced in ifTT) was that the PC 
can be greatly expanded without becoming non-stabilizer. We 
do this by adding Kraus operators that correspond to either 
Clifford operations or Pauli measurements followed by con¬ 
ditional Pauli operations. For the 1-qubit case, the Clifford 
channel (CC), which corresponds to the first expansion, is 
composed of the 24 operators that maintain the symmetry of 
the chiral Clifford octahedron l23l . Another expansion to the 
PC can be obtained by introducing pairs of operators that ef¬ 
fectively produce a measurement in a Pauli basis followed by 
a conditional Pauli operation, such that all states are mapped 
to the same state. We refer to these pairs of operators as 
measurement-induced translations. For each Pauli state, |A), 
these operations can be represented by the following pairs: 

{£ ao = |A)(A|,F; a1 = |A)(A^|} (1) 

The addition of these measurement-induced translations to 
the PC gives rise to the Pauli + measurements channel (PMC), 
while their addition to the CC produces the Clifford + mea¬ 
surements channel (CMC). 


B. Non-stabilizer error channels 


As in ED, we use the efficiently simulable channels to ap¬ 
proximate realistic non-stabilizer channels. In particular, we 
focus on the amplitude damping channel (ADC) and the po¬ 
larization along an axis in the x-y plane of the Bloch sphere 
(Pol^C), shown on equations [2] and [3] respectively: 


Eao = |0)(0| + VT=7|1)<1| 

E A i = Sy\0}(1\ 


( 2 ) 


P04C = 


E X yO — P(p I 

E xy i = y/p^ [cos(» X + sin(0) Y] 


(3) 


models. As we do not want to underestimate the deleterious 
effect of the target channel on quantum information (i.e. we 
want our approximation to be “honest’ ’ ED), we perform the 
minimizations such that the approximate channels constitute 
an upper bound on the error induced on the system. Mathe¬ 
matically, this condition can be enforced in a variety of ways, 
by employing different fidelity or distance measures. From 
the constraints that we have studied, the most lenient one cor¬ 
responds to the average fidelity constraint, in which we en¬ 
force the following condition: 

F av (J, Target) ^ F„(I, Model) (4) 

The average fidelity between a unitary transformation V and 
a quantum channel K is given by: 

FUV , K ) = ^ Y 1 l T r (^)| 2 ( 5 ) 

i 

where N is the dimension of the Hilbert space and { K,} are 
these Kraus operators of the error channel I \. On all the 
approximations that we have performed, the average fidelity 
constraint has always given the same results as if no constraint 
had been applied. 

On the other hand, the most stringent constraint corre¬ 
sponds to the worst trace distance one, in which we enforce 
that for every initial pure state its trace distance to the result¬ 
ing state after the target transformation is not greater than its 
trace distance to the resulting state after the model approxi¬ 
mation: 


E Tr (p, Target(p)) < D Tr (p, Model (p)) (6) 

where the trace distance is calculated using the following ex¬ 
pression: 


D Tl (p,cr) = ^Tr\p-a\ (7) 

This worst trace distance constraint results in approximate 
channels that are honest, in the sense that the deleterious ef¬ 
fect of the target error channel on any pure quantum state will 
never be underestimated, as pointed out by Magesan et al. ua 
and Puzzuoli et al. ED- We could think of an even tighter 
constraint in which we enforce this condition on every initial 
state, pure or mixed. However, if the target and the model 
transformations have different fixed points, this condition is 
impossible to satisfy. We compare results for both constraints 
and label them “a” (average fidelity constraint) and “w” (worst 
trace distance constraint). 


C. Constraints 

The free parameters in the PC and its expansions corre¬ 
spond to the probabilities associated with the Kraus operators. 
Previously, we have obtained stabilizer approximations to the 
realistic non-stabilizer channels by minimizing the Hilbert- 
Schmidt distance ll24l [25l over the parameter space of the 


HI. HONESTY AND ACCURACY AT THE PHYSICAL AND 
LOGICAL LEVELS 

For each approximate channel, we study two properties: 
honesty and accuracy. An approximate channel is honest if 
it does not underestimate the detrimental effect of the target 
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FIG. 1. We study the honesty and accuracy of the approximate channels at the physical level and at the logical level before and after error 
correction. The logical initial state is encoded without errors and then errors are applied. The preparation of the ancillary cat state in the faulty 
EC is error-free as well, as shown in Fig. [2] 


error channel. The accuracy of an approximate channel refers 
to how closely it can mimic the effect of the target channel 
on an initial state. More explicitly, if a target error channel E 
maps a pure state p to E(p) and an approximate channel A 
maps the same state to A(p ), then A is honest if 

D Tt (p,E(p))<D Tr (p,A(p)) (8) 

for every pure state in the initial physical or logical space. The 
accuracy is measured by the average trace distance between 
the resulting states: 

( 'd^(E(p),A(p ))) (9) 

Notice that for both properties, our measure of choice is the 
trace distance. A good approximate channel will be honest (or 
as least dishonest as possible) and as accurate as possible, not 
only at the physical level, but also at the various logical levels. 
We distinguish 4 different scenarios to compare honesty and 
accuracy: (a) the physical (1-qubit) level, (b) the uncorrected 
logical level, (c) the logical level with perfect EC, and (d) the 
logical level with faulty EC, as depicted in Fig. |T] 

For each target non-Clifford error channel, we study 
two kinds of approximations: (a) the Pauli channels (PC), 
which employ only 1-qubit Pauli operators, and (b) the ex¬ 
panded channels or Clifford+measurements channels (CMC), 
which include all the 1-qubit Clifford operators and the 
measurement-induced translations HD. In turn, each kind 
of approximation is performed with the average fidelity con¬ 
straint (“a”) and the worst trace distance constraint (“w”), 
resulting in four approximate channels. Notice that the un¬ 
constrained PC is equivalent to the Pauli Twirled Approxima¬ 
tion Go] ED El, the channel obtained by removing the off- 
diagonal elements from the target channel’s process matrix in 
the Pauli basis |26l . We also analyze the completely isotropic 
Pauli channel or depolarizing channel (DC), the most com¬ 
mon error model used when calculating thresholds. In this 
paper we are comparing single qubit error channels and we 
only use the single qubit depolarizing channel. This channel 
is a version of the PC a where the coefficients corresponding 
to each Pauli matrix are forced to have the same value. This 


error model serves as a reference. The approximations are 
summarized in Table Q] 

TABLE I. Summary of the various target and approximate channels. 


Channel 

Complete name 

Honesty constrained 

ADC 

amplitude damping 

- 

PolC 

polarization along non-Clifford axis 

- 

PC a 

Pauli 

no 

PC w 

Pauli 

yes 

CMC a 

Clifford+measurements 

no 

CMC w 

Clifford+measurements 

yes 

DC 

Depolarizing channel 

no 


IV. CALCULATION OF THE PSEUDO-THRESHOLD 

Our objective with respect to the pseudo-threshold is 
twofold. On the one hand, we want to study how sensitive 
a QECC’s threshold is to the noise model. On the other hand, 
we want to determine if the thresholds obtained with our ex¬ 
panded error models approximate the realistic threshold more 
accurately than the PC. 


A. Procedure to compute the level-1 pseudo-threshold 

Because our target error models are non-stabilizer, we per¬ 
form exact (full density matrix) simulations of QEC circuits 
up to the first level of encoding. We calculate a particular 
QECC’s level-1 pseudo-threshold under a given error channel 
in the following way: 

1. Run the physical circuit: 

(a) Choose an initial 1-qubit pure state, \ip). 
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(b) Apply the selected error channel. 

(c) Compute the fidelity between the initial and final 
states. 

2. Run the logical circuit: 

(a) Encode the initial state using the selected QECC. 

(b) Apply the error channel to each physical qubit. 

(c) Perform EC. 

(d) Compute the fidelity between the initial and final 
logical states. 

We are interested in how much the final state is af¬ 
fected by errors which are uncorrectable by the selected 
QECC. Therefore, for the faultily corrected case, before 
computing the fidelity, we perform one round of perfect 
EC. This has the effect of eliminating correctable er¬ 
rors which happened during or after the faulty EC. The 
process of performing a round of perfect EC and then 
computing the fidelity can also be viewed as computing 
an error-corrected fidelity: 


B. Methods of error correction 

The EC step is performed by measuring the stabilizer gener¬ 
ators and later correcting any detected errors. We distinguish 
between the error-free EC, which results in a code-capacity 
pseudo-threshold, and the faulty EC, which results in the more 
realistic circuit-based pseudo-threshold. The faulty EC is built 
by inserting an error channel after each gate in the original 
circuit. As the Steane code will be the focus of our analy¬ 
sis, consider, for example, the measurement of the stabilizer 
III}.(XXX, as depicted in Fig. [2] The error-free EC step 
would consist of circuits analogous to (a) for each stabilizer 
generator. On the other hand, in the faulty EC regime, each 
stabilizer generator would be measured as shown in (b). Each 
stabilizer measurement is then repeated and the syndrome is 
compared to the one in the previous round. If there is a dis¬ 
agreement between these two, a third round of stabilizer mea¬ 
surements is performed and its syndrome is selected as the 
definitive one. 


V. RESULTS 


Fec(\iPl) , Pl) 


(i/’l I E\ p L PiEi\ip L ) 


1 


( 10 ) 


where is the initial logical state and p l is the fi¬ 
nal logical state, which may not be pure. The set {£)} 
consists of all error operators which the QECC is de¬ 
signed to correct, while { P ,} is the set of projectors to 
the subspaces associated with each error. For the Steane 
[[7,1,3]] code, the set of correctable errors is formed by 
the 64 Pauli operators formed by all possible combina¬ 
tions of X and Z errors acting independently on at most 
one qubit, which includes the Identity operator for the 
case of no errors. 

The set {Ei} includes all (error) operators which the 
QECC is designed to correct. When acting on a given 
state, \4>l), in the codespace, each one of these opera¬ 
tors will transform it into the equivalent state on each 
logical subspace. For example, for the Steane [[7,1,3]] 
code, the set of correctable errors are the 64 Pauli oper¬ 
ators formed by all possible combinations of A' and Z 
errors acting independently on at most one qubit, which 
includes the Identity operator for the case of no errors. 

3. Repeat steps (1) and (2) for various noise strengths to 
obtain fidelities for the physical and logical circuits. 
The threshold is given by the first intersection between 
the two curves. 

4. Repeat this procedure for several initial states to obtain 
an average threshold. For the perfectly corrected case, 
we select 80 initial points uniformly distributed on the 
Bloch sphere. For the faultily corrected one, we select 
20 . 


A. Honesty and accuracy of the approximations 


By construction, the “w” approximations are honest at the 
physical ( 1 -qubit) level, provided that the initial state is pure. 
In our previous work we also determined that when approx¬ 
imating a general non-Clifford channel at the physical level, 
the expanded channels are more accurate than the Pauli. Be¬ 
fore computing the level -1 pseudo-thresholds for different ap¬ 
proximations, we first examine if the honesty of the “w” ap¬ 
proximations and the greater accuracy of the expanded chan¬ 
nels were maintained at the logical level. 


1. Amplitude Damping Channel (ADC) 


For the physical, logical uncorrected, and logical with per¬ 
fect EC levels, we have selected 80 initial states uniformly dis¬ 
tributed over the Bloch sphere surface. For the logical faultily 
corrected level, we have selected 20 points, as the simulations 
involve 3 extra qubits and consequently take an exponentially 
longer time. We have computed the trace distance between 
each one of them and the resulting final state after the ADC 
and its approximations. The average distances are shown in 
the first row of Fig. [3] as a function of the damping strength, 
7 . Likewise, we have computed the trace distance between 
each final state after the ADC and each final state after every 
approximate channel. The average distances are presented in 
the second row of Fig. [3] The behavior in the limit of small 
damping strength (7 —> 0 ) is summarized in Tables [II] and 
[HI] In this limit, it is useful to Taylor-expand the distances 
in terms of the noise strength and compare the coefficients of 
the leading order terms. Expectedly, for the corrected logi- 



data 

qubits 


data 

qubits 


—€ 

3— 

Cl 

rv 


% 

7 




% 

J 

(T) 



•- 



i(|OOM) 

+ 11111 )) 



-( 

He}- 





7 L^J 


TF1, 

._n+i_ 



7 | rj | 

s' H 

S' 

l_EJ Li±J 

fTF 




l “ 1++J 



YXA L+J 

-®- 

HU- 


FIG. 2. Circuits representing the measurement of the operator IIIXXXX in an (a) error-free regime and a (b) faulty regime. In the former 
case, we only need to employ 1 ancillary qubit. Notice that the ancillary qubit starts in the 1+) state and the measurement is performed in the 
X basis. In the faulty EC regime, in order to make the procedure fault tolerant, we employ 4 ancillary qubits initialized in a cat state (27). We 
then measure each ancillary qubit in the X basis and compute their parity to extract the outcome. 


cal cases the linear term is suppressed and the leading order 
is quadratic. For the physical and uncorrected logical cases, 
the leading order is linear. At the logical level with faulty EC, 
simulations were only carried out at low damping strengths 
(7 € [ 10 ~ 5 , 10 -3 ]), which is the pertinent region for the 
pseudo-threshold computation. 

Notice that at the physical level in the first row of Fig. [3j the 
“w” approximations result in curves that are above the target 
ADC by construction, while the “a” approximations produce 
curves below it. This behavior is also present in the small 
noise strength limit, as can be seen by the magnitudes of the 
linear coefficients (Table [Till: PC a < CMC a < ADC < CMC 
w < PC w. Likewise, the accuracies of the CMC approxi¬ 
mations are much better than that of the PC approximations 
(Table |HI]>. In the 7 -+ 0 limit, the CMC approximations are 
w 3 times more accurate. 


At the three logical levels, the “w” approximations are hon¬ 
est for every damping strength. This is true not just on aver¬ 
age, but for every initial state considered. This is an important 
result, as it means that we can safely use the “w” approxi¬ 
mations as a substitute of the ADC when determining codes’ 
thresholds or other error-correcting properties. Remarkably, 
the dishonesty of the PC a is greatly reduced from the phys¬ 
ical to the logical levels in the limit of small 7 . Its error is 
below the honesty cutoff by 36% at the physical level but by 
less than only 2 % for both corrected logical levels and well 
within the deviation in the distance. In contrast, the dishon¬ 
esty of the CMC a is not improved at the logical levels and is 
below the honesty cutoff by 8 — 20 % for all cases. 


The variation of the accuracy from level to level shows an 
interesting behavior. For both levels where the effect of the 
errors is linear (physical and uncorrected logical), in general 
the CMC channels and the “a” approximations are more accu¬ 
rate than the PC channels and the “w” approximations, respec¬ 
tively. This is seen by the magnitudes of the linear coefficients 
(Table III 1 : CMC a < CMC w < PC a < PC w. At the log¬ 


ical level with perfect EC, this intuitively expected behavior 
is seen only for high damping strengths (7 > 0.5) (see Fig. 
|3j. Surprisingly, for lower damping strengths, the most accu¬ 
rate approximation is given by the unconstrained PC, as can 
be observed by the suppresion of the second order terms in the 
accuracy ( Table III 1 . This behavior is particularly pronounced 
at the logical level with perfect EC, where the second order 
terms for the PC a and ADC are practically indistinguishable. 


2. Polarization along a non-Clifford Axis Channel (PoffC) 


We perform an analogous analysis for our second target er¬ 
ror channel: the polarization along a non-Clifford axis on the 
XY plane of the Bloch sphere. We select the axis forming 
an angle (f> = n/8 with respect to the X axis, as this is the 
angle for which the expanded error models perform the worst 
mi- Once again, we have selected 20 initial states for the 
faultily corrected level and 80 points for all other levels. We 
have computed the trace distance between each one of them 
and the resulting final state after the Pol^/gC and its approx¬ 
imations. The average distances are shown in the first row of 
Fig. [4]as a function of the noise strength, p. Likewise, we have 
computed the trace distance between each final state after the 
Pol^/gC and each final state after every approximate channel. 
The average distances are presented in the second row of Fig. 
[4] The behavior in the limit of small noise strength (p -A 0) is 
summarized in Tables IV and[V] As for the ADC, at the phys¬ 
ical and uncorrected logical levels, the leading order is linear. 
At the corrected logical levels, the leading order is quadratic. 


As observed on the ADC, the “w” approximations are hon¬ 
est at every level and for every noise strength. This holds in 
the average case and also for each initial state considered. In¬ 
terestingly, the CMC a becomes honest on average and the PC 
a average distances are indistinguishable from honest. 
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TABLE II. Honesty of the approximate channels to the ADC in the limit of small damping strength. Standard deviations below 10 9 are not 
reported. 


Channel 

Physical 

Logical uncorrected 

Logical perfectly corrected 

Logical faultily corrected 

D Tr /7 

D Tr / 7 

D t V 7 2 

D T V(10 2 7 2 ) 

ADC 

0.55(27) 

3.62 

3.76(96) 

8.0(1.8) 

PC a 

0.347(79) 

3.50 

3.76(96) 

7.8(1.8) 

PC w 

0.81(12) 

8.35 

18.5(3.5) 

37.7(8.0) 

CMC a 

0.50(18) 

4.00 

3.48(45) 

6.3(1.2) 

CMC w 

0.66(24) 

5.33 

6.19(80) 

11.3(2.2) 

DC 

0.333 

3.50 

2.75(36) 

4.95(96) 


TABLE III. Accuracy of the approximate channels to the ADC in the limit of small damping strength. Standard deviations below 10 9 are not 
reported. 


Channel 

Physical 

Logical uncorrected 

Logical perfectly corrected 

Logical faultily corrected 

D Tr /7 

D Tr / 7 

D t V 7 2 

D t 7(10 2 7 2 ) 

PC a 

0.500 

2.41 

7(12) x 10~ 6 

0.123(28) 

PC w 

0.63(26) 

4.94 

14.8(2.6) 

29.8(6.2) 

CMC a 

0.166(60) 

1.35 

1.61(44) 

2.15(74) 

CMC w 

0.194(60) 

1.75 

3.05(94) 

3.7(1.1) 

DC 

0.505(97) 

2.92 

1.68(69) 

3.2(1.2) 


TABLE IV. Honesty of the approximate channels to the Pol„./gC in the limit of small noise strength. Standard deviations below 10 9 are not 
reported. 


Channel 

Physical 

Logical uncorrected 

Logical perfectly corrected 

Logical faultily corrected 

d t 7p 

d t 7p 

d t 7 p 2 

D T 7(10 3 p 2 ) 

Pol^/gC 

0.78(24) 

7.00 

16.2(4.7) 

2.22(72) 

PC a 

0.73(20) 

7.00 

16.2(4.7) 

2.22(72) 

PC w 

0.93(18) 

9.47 

18.5(4.9) 

2.71(77) 

CMC a 

0.77(22) 

7.41 

17.8(5.1) 

2.67(83) 

CMC w 

0.84(23) 

8.17 

20.8(5.9) 

3.12(96) 

DC 

0.667 

7.00 

11.0(1.4) 

1.82(35) 


TABLE V. Accuracy of the approximate channels to the Pol^/gC in the limit of small noise strength. Standard deviations below 10 9 are not 
reported. 


Channel 

Physical 

Logical uncorrected 

Logical perfectly corrected 

Logical faultily corrected 

d t 7p 

D T 7 P 

D Tr /p 2 

D Tr /(10V) 

PC a 

0.278(82) 

2.47 

6.1(7.7) x 10" 6 

5.8(1.6) x 10“ 5 

PC w 

0.36(21) 

4.95 

3.8(1.1) 

6.2(1.6) 

CMC a 

0.108(25) 

1.09 

1.76(41) 

4.7(1.0) 

CMC w 

0.132(43) 

1.37 

4.6(1.1) 

9.2(2.4) 

DC 

0.46(11) 

4.67 

9.2(1.3) 

9.7(1.5) 
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Notice that just like for the ADC, at the physical and uncor¬ 
rected logical levels, the CMC channels and the “a” approxi¬ 
mations are more accurate than their counterparts PC and “w”, 
respectively, as can be seen by the magnitudes of the linear co¬ 
efficients (Table |V): CMC a < CMC w < PC a < PC w. At 
the physical level, and in the p —> 0 limit, the CMC approx¬ 
imations are ss 3 times more accurate than the PC approx¬ 
imations. At the corrected logical levels, the most accurate 
approximation is once again given by the PC a. Surprisingly, 
this behavior holds even up to high noise strengths (p = 0.6). 
In the low noise limit, and at the corrected logical levels, the 
second order terms are practically suppressed. In this limit, 
the PC a is on average more accurate than the CMC channels 
by a factor of 10 5 . 


3. High accuracy of the unconstrained PC 


For both the ADC and the Pol^/gC, the unconstrained PC 
results in approximations that are honest (or almost honest) 
and extremely accurate at the logical corrected levels. In the 
limit of small error, this is very evident by comparing the 
quadratic coefficients of the PC a to the other approxima¬ 
tions (See Tables III and[V]). The high accuracy of the un¬ 
constrained PC in the context of EC has previously been ob¬ 
served. Geller and Zhou found very good agreement between 
the PC a or Pauli twirled approximation and two different real¬ 
istic noise models when correcting a Bell state l22l . Likewise, 
Puzzuoli et al. observed great accuracy of the PC a when cor¬ 
recting a Choi state encoded in the [[5,1,3]] code ED. As 
clearly explained in this Ref. ED, after (perfect) EC, the pro¬ 
cess matrix elements corresponding to Pauli error strings that 
result in different syndromes become zero. Intuitively, we can 
say that the “non-Pauli” advantage of the expanded approxi¬ 
mations at the physical level gets “washed away” after EC. 


imating the pseudo-threshold for every initial state, but result 
in an average that is close to the target’s average. 


TABLE VI. Thresholds for the Steane code under the ADC and its 
Pauli and expanded approximations. ADC/PC a uses ADC at the 
physical level and PC a at the logical level. 


Channel 

Code capacity 

Circuit-based 

(7th) 

RMS 

(7th) x 10 4 

RMS x 10 4 

ADC 

0.18(17) 

- 

4.8(4.2) 

- 

PC a 

0.132(38) 

0.171 

4.8(1.4) 

3.91 

PC w 

0.061(43) 

0.204 

2.36(60) 

4.69 

CMC a 

0.19(17) 

0.0498 

6.4(4.2) 

1.67 

CMC w 

0.15(14) 

0.0644 

4.8(3.1) 

1.12 

DC 

0.162(22) 

0.165 

7.2(1.4) 

4.60 

ADC/PC a 

0.30(37) 

0.255 

4.9(4.2) 

0.101 


TABLE VII. Thresholds for the Steane code under the PoL/gC and 
its Pauli and expanded approximations. PoL/ 8 C/PC a uses PoLygC 
at the physical level and PC a at the logical level. 


Channel 

Code capacity 

Circuit-based 

( Pth ) 

RMS 

(Pth> x 10 4 

RMS x 10 4 

o 

if - 

00 

n 

0.14(24) 

- 

3.5(1.5) 

- 

PC a 

0.086(74) 

0.238 

3.10(26) 

1.53 

PC w 

0.078(16) 

0.237 

3.46(35) 

1.48 

CMC a 

0.11(21) 

0.112 

3.09(85) 

0.816 

CMC w 

0.09(14) 

0.169 

2.91(76) 

0.991 

DC 

0.083(12) 

0.240 

3.92(64) 

1.60 

PoL/gC/PC a 

0.14(25) 

0.0255 

3.5(1.5) 

1.19 x 10“ 3 


B. Level-1 pseudo-thresholds 


We perform the simulations of two different scenarios: (1) 
one with perfect EC, which results in a relatively high code¬ 
capacity pseudo-threshold and (2) one with faulty EC, which 
results in a more realistic circuit-based pseudo-threshold. 
Apart from the average pseudo-threshold, we also calculate 
the root mean square difference (RMS) between the pseudo¬ 
thresholds given by the target non-Clifford channel and the 
ones predicted by each approximate channel: 

rms = V(( P r nei -pr ox ) 2 ) a i) 

The RMS quantifies the accuracy of each approximate 
channel to estimate the pseudo-threshold of the target chan¬ 
nel. We calculate the RMS because comparing only the aver¬ 
age values does not account for any cancellation of errors. A 
certain approximate channel can do a very poor job at approx- 


The results for the ADC and the Pol^/gC are summarized 
in Tables VI and |VII| respectively. For the ADC, the pseudo¬ 
thresholds are expressed in terms of the damping strength, 
7, while for the Pol„./ 8 C, they are expressed in terms of the 
noise strength, p. In both cases, the standard deviation of 
the pseudo-thresholds is included inside parentheses. Notice 
that the code-capacity pseudo-thresholds are about 3 orders of 
magnitude higher than the circuit-based ones. The latter ones 
are on the range expected for the Steane code l28tt . Although 
the code-capacity pseudo-thresholds are unrealistically high, 
they show similar trends with respect to their circuit-based 
counterparts. 


In general, the target non-Clifford channels and their CMC 
approximations result in high standard deviations of the 
pseudo-thresholds around its average values. This implies that 
for these channels the pseudo-thresholds are much more de¬ 
pendent on the initial state than for the PCs. It is interest¬ 
ing that the Pauli channels always result in pseudo-thresholds 
that are lower than the real ones. This trend has also been 
observed by Tomita and Svore on the surface code m and 
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suggests that anistropic Pauli channel approximations to real¬ 
istic noise models are pessimistic. The CMC w approxima¬ 
tions also result in lower pseudothresholds. This is in contrast 
to the isotropic Pauli channel approximation (DC) that yields 
optimistic pseudothresholds. 

The CMCs give more accurate pseudo-threshold estimates 
than the PCs, as can be seen by comparing their RMS values. 
Although we might expect the “a” channels to result in better 
approximations than the “w” channels, in general this is not 
the case. The most important variation is between the CMCs 
and the PCs. In general, however, the “w” channels result in 
lower pseudo-thresholds than the “a” channels, which implies 
that honest approximations at the physical level generate con¬ 
servative estimates of the threshold. Finally, we notice that the 
circuit-based pseudothresholds are quite comparable yielding 
pseudothresholds within a factor of two for all of the error 
models. The DC model representing isotropic depolarizing 
noise yields the least accurate results and the highest thresh¬ 
olds. 

In the previous section, we noticed that the PC a, one of the 
simplest approximations at the physical level, and one that 
is not even honest, results in very accurate and practically 
honest approximations at the corrected logical levels. In the 
context of our level-1 pseudo-threshold estimation, this result 
suggests that we can take a different strategy. Instead of us¬ 
ing the approximate channel at both the physical and logical 
level to calculate the pseudo-threshold, we can use the target 
channel at the physical level and the PC a approximation at 
the logical level. More generally, we can simulate the realistic 
noise model in an exact way whenever it is feasible, and in the 
encoded cases just use the PC a. If we take this approach, we 
obtain more accurate state by state pseudo-thresholds for the 
circuit-based case as seen in Tables IvTl and I VII| 


VI. CONCLUSIONS 

We have studied the feasibility of using approximate error 
channels at the physical level to simulate the performance of 
QEC protocols under the influence of non-stabilizer errors. 
We have selected the Steane [[7,1,3]] code as a model QEC 
protocol and have calculated the honesty and accuracy of the 
Pauli and expanded approximations to realistic non-stabilizer 


errors. We have also computed the code’s pseudothreshold 
under the different error models. 

Similarly to results recently obtained for distance-3 sur¬ 
face codes 0, the PC approximations result in lower pseu¬ 
dothreshold values than the realistic error channels. In con¬ 
trast, the isotropic DC approximation yields higher pseu¬ 
dothresholds than the target channels in the circuit-based 
model. Since most thresholds in the literature use two-qubit 
and one-qubit depolarizing channels, we expect that realistic 
error models with equivalent fidelity will have slightly lower 
pseudothresholds in practice. 

We have also found that physically honest approximations 
compose well: they result in honest approximations at the log¬ 
ical level. Perhaps more interestingly, for both realistic noise 
models analyzed, the dishonesty of the PC a gets greatly re¬ 
duced at the corrected logical levels and its accuracy becomes 
extremely high. This suggests that, if our goal is to estimate 
thresholds, then the best protocol is to model the error as re¬ 
alistically as possible at the physical level and use the PC a at 
the logical level. 

As explained by Puzzuoli et al. ED, single qubit errors can 
be separated into errors that deform the Bloch sphere and er¬ 
rors that preserve the Bloch sphere. Both target error models 
analyzed in our work deform the Bloch sphere where the PC 
a yields effectively honest approximations at the logical level 
for small errors. Puzzuoli et al. found that when the errors are 
unitary, the PC a generally results in a poor approximation. 
In future work, we will compute the level-1 pseudothreshold 
for the Steane and other codes under unitary errors, to de¬ 
termine if Pauli and expanded approximations still result in 
pessimistic pseudothreshold estimates. 
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Physical Logical uncorrected Logical perfectly corrected 



FIG. 3. (color online) Honesty (top) and accuracy (bottom) of the approximate channels to the ADC at various levels. 
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Logical perfectly corrected 




FIG. 4. (color online) Honesty (top) and accuracy (bottom) of the approximate channels to the Pol^/gC at various levels. 













































